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Hot accretion with outflow and thermal conduction 
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Abstract We present self-similar solutions for advec- 
tion -dominated accretion flows with thermal conduc- 
tion in the presence of outflows. Possible effects of out- 
flows on the accretion flow are parametrized and a sat- 
urated form of thermal conduction, as is appropriate 
for the weakly-coUisional regime of interest, is included 
in our model. While the cooling effect of outflows is 
noticeable, thermal conduction provides an extra heat- 
ing source. In comparison to accretion flows without 
winds, we show that the disc rotates faster and be- 
comes cooler because of the angular momentum and 
energy flux which are taking away by the winds. But 
thermal conduction opposes the effects of winds and not 
only decreases the rotational velocity, but increases the 
temperature. However, reduction of the surface density 
and the enhanced accretion velocity are amplified by 
both of the winds and the thermal conduction. We find 
that for stronger outflows, a higher level of saturated 
thermal conduction is needed to significantly modify 
the physical profiles of the accretion flow. 
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1 Introduction 



accretion discs - black hole 



The thin accretion disk model describes flows in which 
the viscous heating of the gas radiates out of the sys- 
tem immediately after generation (Shakura & Sun- 
yaev 1973). However, another kind of accretion has 
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been studied during recent years where radiative en- 
ergy losses are small so that most of the energy is ad- 
vected with the gas. These Adeetion-Dominated Ac- 
cretion Flows (ADAF) occur in two regimes depending 
on the mass accretion rate and the optical depth. At 
very high mass accretion rate, the optical depth be- 
comes very high and the radiation can be trapped in 
the gas. This type of accretion which is known un- 
der the name 'slim accretion disk' has been studied in 
detail by Abramowicz et al. (1988). But when the ac- 
cretion rate is very small and the optical depth is very 
low, we may have another type of accretion (Narayan 
& Yi 1994; Abramowitz et al. 1995; Chen 1995). How- 
ever, numerical simulations of radiatively inefficient ac- 
cretion flows revealed that low viscosity flows are con- 
vectively unstable and convection strongly influences 
the global properties of the flow (e.g., Igumenshchev, 
Abramowicz, & Narayan 2000). Thus, another type of 
accretion flows has been proposed, in which convection 
plays as a dominant mechanism of transporting angular 
momentum and the local released viscous energy (e.g., 
Narayan, Igumenshchev, & Abramowicz 2000). 

This diversity of models tells us that modeling the 
hot accretion flows is a challenging and controversial 
problem. We think, one of the largely neglected physi- 
cal ingredient in this field, is thermal conduction. But a 
few authors tried to study the role of " turbulent" heat 
transport in ADAF- like flows (Honma 1996; Manmoto 
et al. 2000). Since thermal conduction acts to oppose 
the formation of the temperature gradient that causes 
it, one might expect that the temperature and density 
profiles for accretion flows in which thermal conduction 
plays a significant role to appear different compared to 
those flows in which thermal conduction is less effec- 
tive. Just recently, Johnson & Quataert (2007) stud- 
ied the effects of electron thermal conduction on the 
properties of hot accretion flows, under the assumption 
of spherical symmetry. In another interesting analysis. 
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Tanaka & Menou (2006) showed that thermal conduc- 
tion affects the global properties of hot accretion flows 
substantially. They generalized standard ADAF solu- 
tions to include a saturated form of thermal conduction. 
In the second part of their paper, a set of two dimen- 
sional self-similar solutions of ADAFs in the presence 
of thermal conduction has been presented. The role of 
conduction is in providing the extra degree of freedom 
necessary to launch thermal outflows according to their 
2D solutions. 

On the other hand, ADAFs with winds or outflows 
have been studied extensively during recent years, irre- 
spective of possible driven mechanisms of winds. But 
thermal conduction has been neglected in all these 
ADAFs solutions with winds. In advection-dominatcd 
inflow-outflow solutions (ADIOS), it is generally as- 
sumed that the mass flow rate has a power-law depen- 
dence on radius, with the power law index, s, treated 
as a parameter (e.g., Blandford & Begelman 1999; 
Quataert & Narayan 1999; Beckert 2000; Misra & Taam 
2001; Fukue 2004; Xie & Yuan 2008). Beckert (2000) 
presented self-similar solutions for ADAFs with radial 
viscous force in the presence of outflows from the accre- 
tion flow or infall. Turolla & DuUemond (2000) inves- 
tigated how, and to what extent, the inclusion of the 
source of ADAF material affects the Bernoulli num- 
ber and the onset of a wind. In their model, the ac- 
cretion rate decreases with radius (s < 0). Misra & 
Taam (2001) studied the effect of a possible hydrody- 
namical wind on the nature of hot accretion disc so- 
lutions. They showed that their solutions arc locally 
unstable to a new type of instability called wind-driven 
instability, in which the presence of a wind causes the 
disc to be unstable to long-wavelength perturbations of 
the surface density. Kitabatake, Fukue & Matsumoto 
(2002) studied supercritical accretion disc with winds, 
though angular momentum loss of the disc, because of 
the winds, has been neglected. Comparisons with ob- 
servations reveal that the X-ray spectra of such a wind- 
driven self-similar flow, can explain the observed spec- 
tra of black hole candidates in quiescence (Quataert 
& Narayan 1999; Yuan, Markoff & Falcke 2002). Lin, 
Misra & Taam (2001) find that the spectral characteris- 
tics of high-luminosity black hole systems suggest that 
winds may be important for these systems as well. 

Considering extensive works on hot accretion flows 
with winds and significant role of thermal conduction 
in deriving outflows (Tanaka & Menou 2006), we study 
ADAFs with outflows and thermal conduction using 
height-integrated set of equations. A phenomenologi- 
cal way is adopted in which we parameterize the rate 
at which mass, angular momentum and energy are ex- 
tracted by outflow or wind. In the next section, we 



present basic equations of the model. Self-similar solu- 
tions are investigated in section 3. The paper concludes 
with a summary of the results in section 4. 



2 General formulation 

We consider an accretion disc that is axisymmetric and 
geometrically thin, i.e. H/R < 1. Here R and H are, re- 
spectively, the disk radius and the half-thickness. The 
disc is supposed to be turbulent and possesses an ef- 
fective turbulent viscosity. Consider stationary height- 
integrated equations described an accretion flow onto 
a central object of mass M,. The continuity equations 
reads 



dR 



1 dM^ 
27r dR 



0, 



(1) 



where v^i is the accretion velocity (wr < 0) and E — 
2pH is the surface density at a cylindrical radius R. 
Also, p is the midplane density of the disc and the mass 
loss rate by outflow/wind is represented by M„. So, 



]VU{R) = / ATiR'ra^{R')dR', 



(2) 



where mw(i?) is mass loss rate per unit area from each 
disc face. 

Xie & Yuan (2008) derived the height-integrated ac- 
cretion equations including the coupling between the 
inflow and outflow, to investigate the influence of out- 
flow on the dynamics of hot inflow. They showed that 
under reasonable assumptions to the properties of out- 
flow, the main influence of of outflow can be properly 
included by adopting a radius dependent mass accre- 
tion rate. We write the dependence of the accretion 
rate M as follows (e.g., Blandford & Begelman 1999) 



(3) 



where Mq is the mass accretion rate at the outer bound- 
ary Rq. The parameter s describes how the density 
profile and the accretion rate arc modified. In this pa- 
per, typical values of s considered are between s = 
(no winds) and s = 0.3 (moderately strong wind). The 
above prescription for mass accretion rate has been used 
widely for s > (e.g., Quataert & Narayan 1999; Beck- 
ert 2000; Misra & Taam 2001; Fukue 2004) or s < 
(Turolla & Dullemond 2000). Considering equations 
((!]), ^ and we can obtain 



^^A/o( — ) 



(4) 
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The equation of motion in the radial direction is 

.,^|.«,o^-n^)-iA(,a, (5) 

where ft is angular velocity and fix ~ \/ GM^jB? rep- 
resents Keplerian velocity. Also, Cg is sound speed and 
from vertical hydrostatic equilibrium, we have H = 

Similarly, integration over z of the azimuthal equa- 
tion of motion gives (e.g., Knigge 1999) 

RY.v^ — (^Rfl)-—{RvY. — )^^^^, (6) 

where the last term of right hand side represents an- 
gular momentum carried by the outflowing material. 
Here, I — corresponds to a non-rotating wind and 
^ = 1 to outflowing material that carries away the spe- 
ciflc angular momentum it had at the point of ejection 
(Knigge 1999). Also, is a kinematic viscosity coeffi- 
cient and we assume 

V = aCsH, (7) 

where a is a constant less than unity (Shakura & Sun- 
yaev 1973). 

In order to implement thermal conductivity correctly 
it is essential to know whether the mean free path is less 
than (or comparable to) the scale length of the temper- 
ature gradient. For electron mean free path which are 
greater than the scale length of the temperature gradi- 
ent the thermal conductivity is said to 'saturate' and 
the heat flux approaches a limiting value (Cowie & Mc- 
Kee 1977). But when the mean free paths are much less 
than the temperature gradient the heat flux depends on 
the coefflcient of thermal conductivity and the tempera- 
ture gradient. Generally, thermal conduction transfers 
heat so as to oppose the temperature gradient which 
causes the transfer. Tanaka & Mcnou (2006) discussed 
hot accretion likely proceed under weakly-coUisional 
conditions in these systems. Thus, a saturated form 
of " microscopic" thermal conduction is physically well- 
motivated, as we apply in this study. However, one of 
the primary problems for studying the effects of ther- 
mal conduction in plasmas is the unknown value of the 
thermal conductivity. 

Now, we can write the energy equation considering 
energy balance in the system. We assume the gener- 
ated energy due to viscous dissipation and the heat 
conducted into the volume is concerned are balanced 
by the advection cooling and energy loss of outflow. 
Thus, 

^^R del OH- ^'^P ^ a^clR^ dn ^ 



-'^^{RFs)~\w^Am^l{R). (8) 

where the second term on right hand side represents 
energy transfer due to the thermal conduction and 
Fs = 5(/)spCg is the saturated conduction flux on the 
direction of the temperature gradient (Cowie & McKee 
1977). Dimensionless coefficient (ps is less than unity. 
Also, the last term on right hand side of energy equa- 
tion is the energy loss due to wind or outflow (Knigge 
1999). Depending on the energy loss mechanism, di- 
mensionless parameter ry may change. We consider it 
as a free parameter of our model so that larger 77 corre- 
sponds to more energy extraction from the disc because 
of the outflows (Knigge 1999). 

3 Self-similar solutions 

A self-similar solution is not able to describe the global 
behaviour of an accretion flow, because no boundary 
condition has been taken into account. However, as 
long as we are not interested in the behaviour of the 
flow at the boundaries, such a solution describes cor- 
rectly the true solution asymptotically at large radii. 
We assume that each physical quantity can be expressed 
as a power law of the radial distance, i.e. , where the 
power index h' is determined for each physical quantity 
self-consistently. The solutions are 



m)^^o^oi^r-K (9) 

/to 

-R(i?) = --,y^(^)-V2, (11) 
P{R)^c.,Poi^r-K (12) 

(13) 

H{R)^c.RoA, (14) 



itO 

where T,q and Rq provide convenient units with which 
the equations can be written in the non-dimensional 
forms. Note that P{R) is the height-integrated pres- 
sure. By substituting the above self-similar solutions 
into the dynamical equations of the system, we obtain 
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the following system of dimcnsionlcss equations, to be 
solved for lo, wq, Wi, 102 and uj^: 



UJQLU2 ^ m, 



^2 2 1 , /5 



-m - -(s + -)au}3 - si m = 0, 



(15) 
(16) 
(17) 
(18) 



1 
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^^2^3 = -^afujjujs + 5(2 - s) 



(19) 



where m = AIq / {2n RqY.o y/GM^/Ro) is the nondimen- 
sional mass accretion rate. After algebraic manipula- 
tions, we obtain a forth order algebraic equation for 



9«'. l + 2g ,2 4^r5 



l + 2s ,2s 



/(l-2s/2)^ 3 



and the rest of the physical variables are 



2,1- 25^2 
3a 1 + 2s 



/I /5 X 9 1 + 2s 2 4 



1 - 2sP' 



^^2 = ^( 



3q; , 1 + 2s , T 



(20) 



(21) 



(22) 



(23) 



2 n -2sZ2' 
3a 1 + 2s 

We can solve algebraic equation (f20| numerically and 
clearly only real roots which correspond to positive ujf 
are physically acceptable. Without mass outflow and 
thermal conduction, i.e. s = Z = r/ = and 4>s = 0, 
equation (|20p and similarity solutions reduce to the 



standard ADAF solutions (Narayan & Yi 1994). Also, 
in the absence of wind but with thermal conduction, 
equation ([^D|) reduces to equation (17) of Tanaka & 
Menou (2006). But our main algebraic equation in- 
cludes both outflows and thermal conduction. 

Now we can analysis behavior of the solutions in the 
presence of the wind and thermal conduction. Our pri- 
mary goal is to consider the effects of winds and thermal 
conduction via parameters s, I, rj and (ps- First, we can 
summarize typical behavior of the standard ADAF so- 
lutions as follows (Narayan & Yi 1994): (a) The surface 
density increases with the accretion rate, and decreases 
with the viscosity coefficient a; (b) But the radial veloc- 
ity is directly proportional to the viscosity coefficient; 

(c) The gas rotates with sub-Keplerian angular velocity, 
more or less independent of the coefficient a; and finally 

(d) the opening angle of the disc is fixed, independent 
of a and m. Mass outflows and thermal conduction 
may modify these behaviors according to our solutions. 
Equation ((2T|) shows that surface density is directly pro- 
portional to the mass accretion rate. But dependence 
of the surface density on the exponent of accretion rate 
s is determined by solving equation (PU)) . The rota- 
tional and the radial velocities are both independent of 
nondimensional mass accretion rate m. Also, depen- 
dence of the radial velocity on the viscosity coefficient 
is determined by the main algebraic equation. In the 
absence of mass outflows and thermal conduction, the 
opening angle of the disc is independent of the accretion 
rate and the viscosity coefficient and this can be under- 
stood from equation (PO)) by setting s = / = r/ = 0s = O. 
But depending on the angular momentum and energy 
exchanges due to the wind (i.e. I and ry) and the vari- 
able accretion rate (i.e. s) and the thermal conduction 
(i.e. ips) the thickness of the disc may change according 
to the acceptable roots of equation ((20|) . We can see 
that the rotational velocity is sub-Kcplerian according 
to equation ([22|) . 

Figure [1] shows proflles of the physical variables ver- 
sus thermal conduction coefficient 0s for different ac- 
cretion rate exponent, i.e. s = 0, 0.1, 0.2 and 0.3. The 
value of s measures the strength of outflow, and a larger 
s denotes a stronger outflow. The other input parame- 
ters are a = 0.2, 7 = 1.5, / = 1 and I = rj = 1 and each 
curve is labeled by its corresponding s. Recent work by 
Sharma et al. (2006) suggest that the viscosity param- 
eter in a hot accretion flow will be larger than in a stan- 
dard thin disc. If a is much smaller than 0.25, the max- 
imum accretion rate up to which the ADAF solution is 
possible, decreases signiflcantly and the maximum lu- 
minosity of the models becomes much smaller than the 
observed luminosities (Quataert & Narayan 1999). So, 
a = 0.2 is probably not unrealistic for our analysis. 




Fig. 1 Profiles of the physical variables of the accretion disc versus the saturation constant (jis, taking a = 0.2, 7 = 1.5, 
I = 1, Tj = 1 and / = 1. Each curve is labeled by its corresponding exponent s. 




. 2 The same as Figure [T] but s = 0.3 and each curve is labeled by its corresponding coefficient rj. Solid curves are 
1 and dashed lines represent solutions with I — 0. 
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Clearly profiles with s ~ represent no-wind solutions 
(Tanaka & Menou 2006). We can see that ADAFs with 
winds rotate more quickly than those without winds. 
Also, the viscous dissipation is expected to be larger in 
the presence of winds and outflows. Strong-wind mod- 
els have a lower surface density than weak-wind models. 
While solutions without thermal conduction are recov- 
ered at small cps values, we can see significant deviations 
as (j)s increases. Note that for a given set of the input 
parameters, the solutions reach to a non-rotating limit 
at a specific value of (j)s which we denote it by (/fg. We 
can not extend the profiles beyond 0^, because equa- 
tion (j22p gives a negative tol which is clearly unphysi- 
cal. As profiles of Figure [T] show the critical magnitude 
of conduction viscosity 0^ for which the solutions tend 
to non-rotating limit highly depends on the accretion 
rate exponent s. In fact, higher values of s correspond 
to larger (f)"^. 

When there is mass outflow, one can expect lower 
surface density. Figure [T] shows for non-zero s, surface 
density is lower than the standard ADAF solution and 
for stronger outflows, this reduction of the surface den- 
sity is more evident. However, profile of the surface 
density is not affected by (ps when this parameter is 
small. But as conduction coefficient (f>s increases, sur- 
face density decreases. In particular, this reduction due 
to the thermal conduction is more significant for larger 
accretion rate exponent s. In other words, when the 
disc losses more mass because of the winds or outfiows, 
surface density decreases more significantly due to the 
thermal conduction in comparison to the no- wind solu- 
tions. 

Actually, outflows play as a cooling agent and ther- 
mal conduction provides extra heating and there is a 
competition between these physical factors. While ef- 
fects of winds and conduction on the profiles of the sur- 
face density and the accretion velocity are similar, their 
effects on the rotational velocity and the temperature 
oppose each other. Both winds and thermal conduc- 
tion lead to enhanced accretion velocity and reduced 
surface density. Although winds increases rotational 
velocity and decreases the temperature, thermal con- 
duction not only decreases rotational velocity but in- 
creases the temperature too. 

Angular momentum conservation implies (l—2sP) > 
which is trivially valid for non-rotating winds (i.e. 
^ = 0). But for rotating winds, this inequality implies 
s < 1/2. Figure [2] show physical profiles for various r) 
and I, taking s = 0.3, a — 0.2, 7 = 1.5 and / = 1. 
Profiles corresponding to the rotating winds (i.e I = 1) 
are shown by solid curves, but dashed curves show non- 
rotating solutions {I = 0). Each curve is labeled by 
its corresponding coefficient j]. Higher this parameter, 



more energy flux is taking away by winds. Solutions 
with rotating winds are more sensitive to the variations 
of the parameter 77 according to the plots of Figure [21 
As more energy flux is extracted from the disc by winds 
(i.e. higher 77), a lower level of the saturated thermal 
conduction is enough to signiflcantly modify physical 
profiles. In particular, this behavior is more evident for 
solutions with rotating winds. 

4 Discussion and Summary 

Theoretical arguments and observations suggest that 
mass loss via winds may be important in sub - 
Eddington, radiatively inefficient, accretion flows. On 
the other hand, thermal conduction may play a sig- 
nificant role in such systems (Tanaka & Menou 2006). 
Using a simplified model, we included both winds and 
thermal conduction in unified model in order to under- 
stand their possible combined effects on the dynamics 
of the system. Accounting for a variable mass accre- 
tion rate in the flow in proportion to i?" and a saturated 
form of thermal conduction with coefficient (jjs, a set of 
self-similar solutions are presented in this study. We 
have varied s and 0s in our models to judge their sensi- 
tivity to these parameters. The most important finding 
of our analysis is that as more mass, angular momen- 
tum and energy flux are taking away from the disc (i.e. 
stronger winds), possible modiflcations of the physical 
profiles due to the thermal conduction will occur at 
higher level of conduction. 

Although our self-similar solutions are too simplified 
to be used to calculate the emitted spectrum, their typi- 
cal behavior show importance of thermal conduction in 
the presence of winds. For future the global solution 
rather than the self-similar solutions is required. This 
is because most of the radiation of an ADAF comes 
from its innermost region where the self-similar solu- 
tion breaks down. 
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